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•'-H . The light ainphfication by finite active media is used extensively 

J>-„ in modern optics applications. In this paper the light amplification 

"^ I and scattering by cluster of small particles is studied analytically and 

numerically with the help of the local perturbation method and phe- 

nomenological laser theory. It is shown that light amplification is pos- 

^ ■ sible even for one small particle, and that the amplification is more 

QQ ' profound when the light frequency nears the frequency of the cluster's 

^N . morphological resonance. Theoretical discussions are supplemented by 

. . numerical results for scattering by clusters which particles positioned 

^^ I at ordered and at slightly disordered positions. 



1 Introduction 

Light amplification by stimulated emission of radiation (laser) is well known 
phenomenon routinely employed in solids-state, gas, and dye lasers PP-[3]. 
With the development of micro and nano technology the new kind of lasers, 
so-called photonic crystal lasers, emerged |1]-[Z]. The common peculiarity of 
the photonic crystal lasers is that they are made of finite number of parti- 
cles or cells ordered at least in one dimension. These lasers are extremely 
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versatile: they can be made from active host medium filled by passive scat- 
tering particles, or from active scatterers immersed into passive host medium, 
or both. Moreover, the particles in such lasers can move or be static. From 
practical point of view, the photonic crystal lasers and amplifiers can be used 
as light sources to compensate the optical losses in metamaterials [8]-[T0]. 

To the author's knowledge, the light generation by a scattering medium 
with negative resonance absorption was initially studied in work [TT] where 
it was shown that the lasing is possible in such medium. The scattering by 
individual active spheres and cylinders was studied analytically and numeri- 
cally in a number of works (see, for example [12]- [15] and references therein). 
At the same time, the discussion of the scattering by a cluster of active par- 
ticles is somewhat limited to very small clusters made of few particles |16] . 
Scattering by many active particles occurs in so-called random lasers ([17]- 
[22]), while analytical predictions are difficult to made for such systems due 
to extremely large number of particles. 

Recently, the scalar wave scattering by dispersive particles was studied 
by using the local perturbation method (LPM) in the work [23]. In reality, 
however, the vector wave scattering occurs. 

To the author's knowledge, the cluster amplifier made of small particles 
was not studied analytically and numerically in the literature before. 

In this paper the light amplification and scattering by active particles with 
the size smaller than the incident wavelength is studied with the help of the 
LPM [23]- [27] and the phenomenological laser theory [I]- [3]. As an example, 
two important cases are studied numerically: scattering by ordered cluster 
and by weakly disordered cluster. It was shown that the light amplification 
is possible in both cases, however it is severely affected by the size of the 
particles, by the concentration of the doped active atoms, by interaction 
between the particles, and by morphological resonances. 

2 The LPM formalism 

The formalism used in this section is described in a number of papers (|24j- 
|25j ) and it will be briefly presented here for convenience and consistency. 
Consider the cluster positioned at the origin of the coordinates and made 
of N identical active particles which characteristic size L is small compared 
to the incident wavelength A. The electric field E propagating in the host 
medium filled with the particles is described by the following equation [26] 



(A - V 8) V + A;2) E(r) + - V /(r - r„)(£sc,n - 5,)E(rO = S(r), (1) 

where r and r„ are the radius vectors of the observer and the n-th particle 
respectively, and 

Here A and V are the Laplacian and nabla operators, ® defines tensor 
product, fc is a wave number in the host medium (a; is the angular frequency 
and c is the speed of light in vacuum), esc,n and Eh are the relative (in 
respect to vacuum) permittivities of the n-th particle and the host medium 
respectively, / is the function describing the shape of the scatterers, Vn is the 
volume of the n-th particle, and S is the source of the field. The permittivity 
Escn of the active particles, in principle, is complex function depending on the 
electric field E(r„) inside the particle, the frequency u, and other parameters. 
We will discuss this topic in greater detail in the next section. 

It should be noted, that the equation ([1]) is an approximate one and it is 
correct only when the small scatterers [kL <^ 1) are considered. 

The solution of the equation ([T]) can be expressed in the form 

E(r) = E,„(r) + E,,(r), (3) 

where the scattered field E.r is 
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E,,(r) = £ ("/ + ^-^ ) > ;(e,,,„ - e,)E(rO$„(r), (4) 



and 



Here / is the 3x3 unitary tensor in polarization space and r„ is the radius 
vector of the n-th particle. The function / is the Fourier transform of the 
function /. The incident field Ej„ is created by the source S in the host 
medium (more information can be found in j28]). 



The formula (jl]) is rather general one and it describes the field scattered 
by the cluster made of small particles of arbitrary form. The resonance 
properties and interference between the scatterers are taken into account by 
the fields E(r„) inside the particles. The fields E(r„) should be found by 
solving the system of 3A^ linear equations obtained by substituting r = r„ 
into (^. 

The scattered field (jl]) can be simplified even further when the distance 
Rn between the observer and an n-th scatterer is large, i.e. when i?„ ^ L. 
In this case the integrals (|5]) can be calculated explicitly and the scattered 
field (jl]) can be presented in the following form 

^ ^ n=l 

where 

Rn=\r- r„,| > L. (7) 

Here Rn is the distance between the observation point r and the radius vector 
r„ of the n-th scatterer, V is the scatterer's volume. 

In many practical cases the distance between the cluster and the observer 
is much larger than the size of the cluster (|r| 3> max(|r„|)) and the incident 
wavelength (fc |r| ^ 1). In this case the field (j6]) can be simplified and it can 
be rewritten in the following form 



E,e(r) 
where 
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(T-\(S)\^ J2i^sc,n - £hMrn)e~ 



n=l 



1 = r/r, r = |r| ^ max(|r„|), kr ^ 1. (9) 

3 The permittivity of the active particles: steady 
state solution 

In this section we will study the permittivity 6sc,n of the active particles which 
characteristic size is much smaller than the incident wavelength [kL -C 1). 
We note that the permittivity esc,n in formula (J8]) can be complex number 



with negative or positive imaginary part, and in this case one can study wave 
scattering with gain or loss in active media [12]- IE]- The problem with such 
approach is two- fold. First, the value of the imaginary part is not related 
to the properties of the actual medium, and second, the permittivity is the 
same for all particles, that is not true for real systems. As it was suggested 
in [15], such approach is valid for quantitative estimations of lasers before 
threshold. For more accurate investigation one should use rigorous methods 
taking into account atomic transitions and pump dissipation. 

It is important to acknowledge that when the particles are small the 
number of active atoms in upper state is constant within the particle, while 
this number can be different for other particles. The particles are active 
due to the doped active atoms with the density M. We assume that all 
the scatterers have the same density of the active atoms. It is not limiting 
assumption, and it is very convenient one. We also assume that the amplifier 
is at steady state, so that the densities of the active atoms in the upper and 
low states are time independent. 

We approximate the active atoms as two-level systems exited by the op- 
tical pump with the frequency Up and relaxing with the wide range of the 
frequencies ujj. In this case we can present the permittivity esc,n of the par- 
ticles in the following form 

£sc,n = ^sc,n "•" ^sc,nV-ri^ ^j^ l-'-Uj 

where £°^ „ is the permittivity of the n-th scatterer without the active atoms 
and e'g^ „ is the permittivity of the n-th particle due to the presence of the 
active atoms. The latter can be expressed in the following form [IJ 






^4^[Mi(r„)a,(u;) - Mu{v^)a,{uj)] , (11) 

where the emission and the absorption cross sections of the active medium 
respectively are 



and the density of the active atoms is 



M = Muivn) + MLirn). (14) 

Here Mu and Ml are the densities of the active atoms in the upper and 
low states respectively, and M is the total density of the active atoms. The 
frequencies Qa and Qe are the resonance frequencies for absorption and emis- 
sion respectively, and a is the oscillator strength. The frequencies 7^ and 
7g are the dipole relaxation frequencies for the absorption and the emission 
respectively (typical values are about 10^^ Hz [29]), and e and m are the 
electron's charge and mass respectively. 

It should be noted that we neglected by the real part of the permittivity 
e'g^^ because it is small compared to the optical contrast e^^^ — Sh, especially 
near the resonance. We note that the formulae (IT2|) were calculated assuming 

that U ~ Qe, ^a- 

The importance of the formula (ITT]) is that it allows us to use experimen- 
tally measured emission and absorption cross sections (cXe and 0"^) when the 
density Mu is known. 

When the pulse duration exceeds the dipole relaxation time (typically 
10^^'^ s), the density Mu of the upper level atoms can be calculated by using 
the rate equation approximation in which the dopants respond so fast that 
the induced polarization follows the optical field adiabatically [2S]. We use 
the following rate equation [31] 



dMuirn) ^J2lj{rn,u) [ML(r„)a,, - M^(r„)ae,] - Mc7(r„)/r, (15) 



dt 
where 



aaj = (Taiojj), aej = (Teiu;j), (16) 

lA^n,u;,) ^ ^\E{r^,u,)fAu\ (17) 

Here aaj and aej are the absorption and the emission cross sections at the 
frequency Uj respectively, and Ij is the flux of photons at the frequency Uj , 
and the |...| brackets denote the absolute value. The parameter r is the 
relaxation time of the exited atom (typically 10"'^ — 10~^ s [29j), h is the 
reduced Planck constant, and Au = ouj+i — ujj is the frequency bin. 
Since we consider only steady state solutions when 
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f?^%M=o. (18) 

at 

the solution of the rate equation (TTB1) is 

and the formula flTTl) for the permittivity e:^^ „ has the form 

cM r—aa{uj)/T + Y,- Ijir^, ujj) [aa{uj)(Tej - ae{uj)aaj] 



' (r- ,.A = o- /cO 



e,,„{r„,uj = I X e 



(20) 

The formula (l20l) is the main result of this section and it suggests that the 
permittivity of the active scatterer is a complex function of the intensities of 
the fields E(r^,co') inside the particles, frequency u, position of the particles 
r„, and absorption and emission cross sections aaj and aej- When the light 
intensity changes (due to increased reflection from the scatterer's boundaries 
or due to decreased amount of the scattered light from all other particles, for 
example), it will affect the permittivity e'^^^ (I2U|) . We note that the photon 
fluxes Ijijn) should be found by solving the system of nonlinear equations 
with respect to the fields E(r^,a;j) inside the particles. 

Since we consider the scattering by small particles, it is natural to use 
perturbation theory where the Ij is small compared to the pump Ip. When 
the signal is so weak that 

1/r + /p(r„, Up){(Tap + (Tep) > JZ ^■'■(^"' ^i)^^'^o + ^ej), (21) 

O'ap = (^a{.^p), CTep = (^e{^p), (22) 

the permittivity (l20l) can be expressed in the approximate form 



^^^•'^^'""''^^ ' u; V^-''^ l/r + Ipir^,Up)iaap + a,p) 

(23) 
For some estimations it can be sufficient to use the permittivity fl23p . while 
for rigorous numerical calculations one can apply general formula f l20|) where 
fluxes Is are found by using the method of successive approximations. 



When emission and absorption spectra are very distinct and separated 
such that (Jep = (Jas = 0, we can simphfy general formula fl^U]) for the pump 
and the signal respectively 






oOp V ^^''^ l/r + Jp(r„)a, 



ap 



(24) 



,' rr„,.,) = -,£j£./;r ^^'■•■>""°- 



Wo 



l/r + /p(r„)(Tap' 



The important feature of the formulae i^^ is the sign flip: for pump it 
is positive (pump is absorbed) and for the signal it is negative (signal is 
amplified) . 

4 Intensity of the scattered field and the am- 
plification conditions 

We define the intensity of the scattered field as Isc = \Esc\ , and by using the 
formula ([8]) we can present the intensity in the following form 



where 



'""^ IGn'^e'lr'^ 



N 

(T~ I (g) l) ^(e.c,n - eh)E{rn)e- 



2 

, (25) 



1 = r/r, r = |r| ^ max(|r„|), kr ^ 1, (26) 

and the permittivity esc,n is described by the expression (!23|) or by (120|) . 

The expression (12^ suggests that the light amplification (related to the 
imaginary part of the permittivity esc,n) is due to step wise amplification, 
and it is coded in the fields E(r„). Below we consider the fields E(r„) in 
grater detail. 

4.1 The light ampUfication by small sphere 

Consider the light amplification by one small sphere. In this case the intensity 
of the scattered field is described by the expression fl25l) where the field E(ri) 
has the following form [30] 
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E(ri,cj) 



D(uj) 



and the denominator D is 



(27) 



D(cu) = 1 + ^^'''' ^"^ (l - L'^k^ - t-L^e 
^ ' Sen \ 3 

The resonance frequency Ur is found from the following equation 



(2^ 



ReD{uJr) =0, 
and the resonance width ^ is defined as 

Ira D(u)r) 



(29) 



aRcD(a;) 



did 



(30) 



In accordance with formulae (129!) and (l30l) the resonance width and frequency 
respectively are 



e 



and 



OJr 



3c 



^Vs 



+ 



3c^ Im(£:^c,i) 



2a;^L2 Re (e,c,i - ^h) 



(31) 



1 + 



Re(esc,i — ^h] 



1/2 



(32) 



Ly/Re(£:^c,i - Sh) V ' 3e,, 

When the permittivity e^ci of the particle is real, the resonance width is 
defined by the first term in Eq. (I3ip . When the imaginary part of the 
permittivity esc,i is taken into account and it is negative, the resonance width 
can be slightly decreased. This decrease corresponds to effective gain of the 
field scattered by the particle at the resonance frequency. 

The formulae fl^ and fITT]) suggest that when Nua^ < N^aa, the reso- 
nance width increases (with respect to the one in the passive medium) and 
it decreases when Njjae > Niaa- 
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Figure 1: The normalized intensity of the scattered field ^ versus normalized 
frequency kL for the spherical cluster made of small cubes with different 
density of the active atoms M. The permittivity of the particles and the 
host medium is esc,n = 2.53 and Eh = I respectively, the characteristic size of 
the cubes is L = 25 nm, and the total number of the particles is A^ = 123. 
The point source is positioned at r^ = {1,0,0} and the observation point is 
r = r,. 



5 Two numerical examples: light amplifica- 
tion by ordered and weakly disordered spher- 
ical cluster 

In this section the light amplification and scattering by clusters made of 
small active cubes is studied numerically. The normalized intensity 3fJ of the 
scattered field is calculated. The intensity is defined as 

5R(r, u) = Isc{r, uj)/linc{0, u). (33) 

It is assumed that the clusters are made of cubes doped with Yb^^ (active 
material). The absorption and emission cross sections are taken from [32] and 
the other parameters are from [31]. The incident field is generated by the 
point source described by the following formula 
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■run 1, rundom 
run 2, rundom 
■not rundom 
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Figure 2: The normalized intensity of the scattered field ^ versus normalized 
frequency kL for slightly disordered spherical cluster made of small cubes. 
The results of two runs are presented, and the scattering by the ordered 
cluster is also shown for comparison. The density of the active atoms is 
M 
e. 



10 m~ . The permittivity of the particles and the host medium is 
-sc,n = 2.788 and Eh = I respectively, the characteristic size of the cubes 
is L = 25 nm, and the total number of the particles is about A^ = 123. 
The point source is positioned at r^ = {1,0,0} and the observation point is 



E. 



r,u) 



gifc(r-rs) 

Eo^ :, A;|r-r,|>l, (34) 

where the field Eg is polarized along z direction. The source is positioned 
at r^ = {1,0,0} and the center of the cluster is positioned at the origin of 
coordinates r = {0,0,0}. The pump wavelength is 911 nm. At this specific 
wavelength the field Eq is artificially increased by several orders of magnitude. 
Finally, the perturbation theory is used to calculate the density of the upper 
level atoms. 
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5.1 The amplification and scattering by ordered spher- 
ical cluster 

Consider the scattering by the spherical cluster made of the small active 
cubes organized into simple cubic lattice. The size of the cubes in the cluster 
is L = 25 nm, the period is d = 2L (continuous medium), the radius of 
the cluster is 3d. The permittivity of the scatterers is esc,n = 2.53 and the 
permittivity of the host medium is £/i = 1. Such combination of the cluster's 
permittivity and dimensions creates optical resonance of the cluster near 
kL = 0.162 (A = 970 nm). 

Th figure [T] shows the normalized intensity 3fJ of the scattered field for 
the cluster made of the material with density M = 10^^ m~^ (solid red line), 
M = 5* 10^^ m~^ (dashed green line), M = 10^^ m~^ (dash-dotted blue line), 
and M = m~^ (solid black line) for comparison. 

The figure shows that light is profoundly amplified near the frequency 
kL = 0.162, and in addition, the light is weakly amplified at the frequencies 
near kL = 0.158. This feature is probably related to the emission spectrum 
of active material {Yb^~^) which has two crests: very narrow one near 980 nm 
and broad one near 1030 nm. The figure also suggests that at high doping, 
the light amplification increases several orders of magnitude while the doping 
increases only few times. 

We realize that the size of the cluster is too small to made lasing at the 
conventional values of M (~ 10^^ m~'^) and that is why we increased the 
number of the active atoms well above a realistic limit (~ 10^^ m~'^). 

5.2 The amplification and scattering by spherical clus- 
ter with weak positional disorder 

In this subsection we consider the light amplification and scattering by the 
cluster which particles are randomly positioned near predefined positions. 
The predefined positions are the nodes of the cubic lattice with the period 
oi d = 2.2L, and the particles are positioned not further than O.IL from the 
nodes to avoid a collision. We note that the distance O.IL is actually 2.5 nm, 
that is much less than the size of the particle L, and that is why we call this 
cluster weakly disordered one. 

The results of the calculations are presented on the Figure [2j The figure 
shows two typical results for disordered cluster (dashed and dash dotted 
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lines) and one result for ordered cluster (solid line) for comparison. 

The figure suggests that despite weak positional disorder the random 
positioning significantly influences the scattering and amplification by the 
cluster. In our particular case, two scattering peaks compete with each other: 
one near kL = 0.161 (A = 976 nm) and another is near kL = 0.154 (A = 1020 
nm). 

The random amplifier differs from the nonrandom one in a number of 
ways. The first difference is the absence of well defined boundaries, which 
in turn, govern the morphological resonances. Thus, the amplification (or 
lasing) can be at several frequencies simultaneously. The second difference 
is the random structure inside the cluster, affecting the interaction between 
the particles and the total gain as the result. 

6 Conclusions 

The amplification and scattering by the cluster made of the small active 
particles have been studied analytically and numerically. It has been shown 
that such cluster can amplify light when the concentration of the active atoms 
is sufficiently high. It has been shown that the amplification occurs near the 
morphological resonances which are governed by the shape of the cluster and 
its effective optical contrast. 
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